Simulation of bipartite qudit correlations 
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We present a protocol to simulate the quantum correlations of an arbitrary bipartite state, when 
the parties perform a measurement according to two traceless binary observables. We show that 
log(d) bits of classical communication is enough on average, where d is the dimension of both systems. 
To obtain this result, we use the sampling approach for simulating the quantum correlations. We 
discuss how to use this method in the case of qudits. 



I. INTRODUCTION 



o 
o 

In 1964 John Bell showed that the correlations exhibited by the EPR gendanken experiment P, Q could not be 
reproduced by a so-called local hidden variable model, that is, a model where the parties share an infinite amount 
of locally created hidden variables. This nonlocal aspect is one of the strangest properties of quantum physics, and 
OO . understanding this notion remains an important problem. Recently, quantum information processing has provided 
a new point of view to understand quantum nonlocality. In particular, the framework of communication complexity 
has provided tools to study nonlocality. For example, two parties, whom we call Alice and Bob, cannot reproduce 
quantum correlations if they share only hidden random variables (shared randomness), but in some cases, if they are 
allowed to use some additional resources, it becomes possible for them to reproduce the quantum correlations. It is 
precisely this amount of additional resources which we consider here; they allow us to quantify quantum nonlocality. 

The problem of reproducing the statistics of projective measurements on the singlet has been widely studied, with 
communication as the additional resource. In 1992, Maudlin |j| presented a protocol in the case of mesurements 
in the real plane and proved an average-case communication upper bound of 1.17 bits, and independently in 1999, 
Brassard, Cleve and Tapp gave a protocol, together with a worst-case communication upper bound of 8 bits, for 
' arbitrary projective measurements Q. In 2000, Steiner, independently of Maudlin, gave a protocol for projective 
Oh, measurements in the real plane with an average-case upper bound on communication of 1.48 bits, and Cerf, Gisin 
' and Massar [f| proved that for an arbitrary projective measurements, 1.19 bits of communication sufficed on average. 
Recently, in 2003, Toner and Bacon have shown that one bit of communication is always enough to reproduce the 
quantum correlations for arbitrary projective measurements on the singlet state 

Some other resources have been used to simulate quantum correlations resulting from projective measurements on 
the singlet state. These include post-selection 0, and nonlocal boxes In 2005, we have shown that simulating these 
quantum correlations could be reduced to a sampling problem, from which we derived many of the above-mentioned 
, protocols, in a unified framework 0- 

Nevertheless, these results address the simplest scenario, that is, simulating the correlations resulting from mea- 
surements on the singlet state (mostly for projective measurements, with a few extensions to POVMs). There are few 
results about non-maximally entangled pairs, multiparty states, higher dimensional states, or more general measure- 
ments. One significant result in this direction is a protocol from Massar et al able to reproduce the correlations of 
arbitrary measurements on an y en tangled pair of d-dimensional states (qudits) using 0(d\ogd) bits of communication 
but no local hidden variables |1C|. 

In this paper, we use the sampling approach developed in and generalize it to the case of a bipartite pair of 
arbitrary-dimension states (qudits). We study the case where the parties make a restricted type of measurement with 
only two opposite outcomes {1, —1}, that we call traceless binary observable, or TBO. 

Furthermore we impose no constraint on the bipartite (pure) state whose correlations the parties wish to simulate; 
it could be maximally entangled or non-maximally entangled. For an arbitrary bipartite qudit pair, we show that 
log(cZ) bits of communication on average are enough to simulate the joint correlations of the outcomes (where the 
joint correlation is defined as the expectation value of the product of Alice's and Bob's outcome). In the special case 
of maximally entangled qudit pairs, our protocol also reproduces the marginal probabilities, and therefore the full 
probability distribution. 

We will begin by describing the quantum correlations in arbitrary dimensions that we want to simulate classically. 
Then, using the sampling approach, we will present a generalization of the local biased hidden variable model for 
arbitrary dimensions, and present a classical protocol which uses log(c?) bits of communication to simulate the joint 
quantum correlations of an arbitrary bipartite qudit pair. 
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II. THE QUANTUM CORRELATIONS 

In this section we describe the system that we want to simulate classically using some communication. Two parties, 
Alice and Bob, share an arbitrary bipartite qudit pair. They each perform a measurement on their part, where the 
measurements are restricted to what we call traceless binary observables, described below. 

We describe the measurements using observables instead of measurement operators. We restrict the measurements 
to be such that only two outcomes are possible, and these outcomes are equally likely when the measurement is 
applied to a maximally mixed state. 

Definition 1 (Traceless Binary Observable). An observable A is called a Traceless Binary Observable (TBO) if 

• the observable is traceless (i.e. Tr(A) = ) 

• the outputs of the measurement are two opposite values, more precisely A 2 = 11. 

We describe the bipartite quantum correlations on qudits obtained when a TBO is applied to each part of the state. 

Definition 2 (Qudits TBO experiment). Two parties, Alice and Bob, share an arbitrary bipartite qudit pair, 
|V>) € H d ®H d . Alice and Bob measure their part of the state according to their input, describing a TBO pair A,B. 
They then obtain measurement outcomes As {1,-1} and Be {1,-1} respectively. 

We will use the following notation throughout the paper. Let Oa denote the space of TBOs over Tt d . We use §„ to 
denote the unit hypersphere in R Tl+1 . (For example, §2 is the unit sphere in R 3 .) We will also use S n to denote the 
surface area of S n . 

Tsirelson showed that there is a function that maps TBOs over the Hilbert space (matrix operators) to points 
on the surface of a hypersphere, as follows \T\\. This formulation of Tsirelson's theorem was pointed out to us by 
B. Toner ^1 an d appears in this form in |13(. 



Theorem 3 (Tsirelson). For any d > 0, and \ip) E 7i ® H d , there is a function v : Od — > §2d 2 -i such that the 
following holds. If A and B is a TBO pair over TL d , and A, B are the outcomes of measuring according to A and 
B, then 

E(AB\A ® B) = v{A) ■ v(B). 

In the remainder of the paper, we use this theorem implicitly, and use the notation a = v{A), and b = v{B), to 
denote Alice's and Bob's inputs (or measurement), respectively. Furthermore, we call E{AB\A®B) the joint quantum 
correlations. 

III. THE CLASSICAL PROTOCOL 

We present a protocol to simulate bipartite qudit joint quantum correlations as defined in the previous section, 
where the state \ip) is an arbitrary state in Ti d (8 Ti d . First, we generalize the sampling method introduced in 

A. Local hidden biased variable model 

As in @, we consider a model where Alice and Bob share random variables that can depend on Alice's and/or 
Bob's input, which we call a local biased random variable model. 

We generalize the sampling theorem in to arbitrary bipartite qudit states, with TBO measurements. We start 
with a technical lemma to compute the normalization factor of the biased distribution. 

Lemma 4. /„ b- X s dX s = |<5„_i. 
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Proof. Since dX s — d0 n sin(0 n -i)d0 n -i sin 2 (0 n - 2 )d0 n - 2 ...sin n (0 2 )d0 2 sin n 1 (0i)d0i 



b-X. 



d0 n mn(0 n _i)d0 n _i / sin\0 n „ 2 )d0 n _ 2 ... sin n - 2 (0 2 )d0 2 sin n -\0 1 )\cos(0 1 )\ d0x 

iq Jo Jo Jo Jo 

S n - X 2 / sin n - 1 (0 1 )cos(0 1 ) d0i = 2 S n ^ 1 [sin n {0 1 )/n)l' 2 
Jo 

2 s 

— <->ra-l- 



□ 



We will write R n — ^<S rl _i to simplify notation. 



Theorem 5 (Generalized sampling theorem). Let a and b € S n be Alice's and Bob's inputs. If Alice and Bob 
share a random variable X s € §„ distributed according to a biased distribution with probability density 



p(X s \ab) = p s (X s ) 



a ■ A, 



Rn 



then they can simulate the joint correlations 



with marginal expectations 



E(AB\ab) =a-b, 



E(A\ab) = E(B\ab) = 0. 

This says that in this model, Alice and Bob can simulate these correlations without any further resource, that is, 
simulating the bipartite two output joint quantum correlations reduces to distributed sampling from the distribution 
PS- 
Proof Consider the protocol where Alice and Bob set their respective outputs to A(aX s ) = sgn(a • A s ) and B(bX s ) = 
sgn(6 • A s ), where sgn(x) = 1 for x > and sgn(x) = — 1 for x < (x £ M). Then the joint expectation E(AB\ab) is 
given by 



E(AB\ab) 



Pe{X s ) A{a,X s )B{b,Xs) dX s 



a ■ A< 



sgn(a • A s ) sgn(6 • A s ) dX s 



R 

— I (a • X s ) sgn(6 • A s ) dX s 
H 

1 



■n JS„ 



= A s sgn(6 • A s ) dA s ). 

Rn Js n 

Observe that the final integral is invariant by rotation around 6, so it must be the case that 

A s sgn(6 ■ A s ) dX s = cb, 
with c a real constant. 

Multiplying Equation ^ by b on either side to compute the constant, we obtain 

b ■ A s sgn(£> ■ A s ) dX s = c (6 • b) = c. 



(1) 
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By Lemma c = R n , therefore, 



E(AB\ab) = — a-b = a-b. 

Rn 



Finally, we compute the marginal distributions. It is easy to see that 



E(A\ab) 



ps(Xs) A(a,X a ) dX s 



R 



n JS 



'• A< 



sgn(a • A s ) d\ s 



1 

Rn 



(a ■ As) dX s = 0. 



For the second marginal, 



E(B\ab)= ps(As) B(b,X s ) dX s = 



Rn 







L 


d ■ a s 







sgn(6 • A s ) dX s 



With S+, §_ the half-spheres with respect to a, and with A + G § + and A_ £ S_, 



E(B\ab) = 



1 

Rn 



a ■ A+ sgn(& • A+) dX + — 



1 

Rn 



a ■ A_ sgn(6 • A_) dX- 



We make a variable substitution A_ = — A + and we obtain 



R 



E(B\ab) = — a ■ A+ sgn(6 • A+) dA+ - — / a- A+ sgn(6 • A+) rfA+ = 



n JS + 



R 



n JS + 



□ 



Note that we obtain E(A\ab) — E(B\ab) = 0. In the case of measurements on maximally entangled states, it is also 
the case that the marginal expectations are zero. This is because the reduced states of Alice and Bob are maximally 
mixed and, for traceless binary observables, the marginal distributions are uniform, so that E(A\ab) = E(B\ab) = 0. 
However, for arbitrary states, our method will reproduce the joint quantum correlations E(AB\ab), but not necessarily 
the marginal distributions. 



B. Sampling the biased distribution: the rejection method 

It now remains to show how Alice and Bob can obtain a shared sample X s G §„ distributed according to the biased 
distribution Pa(X s ) = a ■ X s / R ni with help of shared A uniformly distributed on §„, and with the additional help 
of communication. Using the same idea as |9L ll4llT5j |. Alice uses the rejection method to perform the sampling. 

Theorem 6. There is a local hidden variable protocol that simulates the joint correlations of TBO measurements on 
a bipartite pair of d- dimensional states using log(c£) + 0(1) bits of communication on average. 

As noted above, for arbitrary states, we reproduce the joint correlations E(AB\ab), and in the special case of 
maximally entangled states, we reproduce the full joint distribution exactly. 

Proof. By Theorem^ it suffices to give a protocol to sample the distribution pg on §„ for n = 2d 2 — 1 in a distributed 
fashion. We show that this can be achieved with log(ci) communication on average. We obtain a sample by applying 
the rejection method [16|, using unbiased (uniform) shared random variables. 

Let U(X) be the uniform probability density function on S n , that is, U(X) = 1/S n . The protocol is as follows. 

1. Alice obtains a uniform sample A ~ U 

2. She computes |a-A|, and accepts A with the corresponding probability. If she accepts, she sends Bob the iteration 
at which this occurred. 
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3. If she rejects, she starts over with a new sample A. 

We compute the probability that Alice accepts at a given iteration (let us call this event "ok" ) , on average over the 
choice of A. 

p(ok) = f p(ok|A)p(A) dX. 
Js„ 

Since A is uniformly distributed on S n , we have p(X) = U(X) = X/S n . Moreover, we accept a given A with probability 
p(ok|A) = \a- A|, so that 

P(ok) = 4" / \a-X\dX = ^, (2) 

where we have used Lemma 0] 

We may now compute the distribution of accepted A's as follows 

/?, , \ P(A) p(ok|A) 
p(A| ° k) = p(ok) 

_ |« - A| 



which corresponds to as required. This proves that the protocol achieves its goal. It remains to show that it 
requires at most 0(log(ci)) bits of communication on average. 

The message sent in the protocol is the iteration i at which Alice accepts the current uniform sample. The 
distribution of the messages behaves according to a Poisson distribution P p , that is, i is sent with probability 

p p {i) = {i- P y- 1 P , 

where p = p(ok) in our case. To compute the number of bits sent on average, it suffices to compute the entropy of 
this distribution, which we do below. From Equation [3 and Lemma we have 

2 Sn-l 

p(ok) = ^s7- 

The surface area S n of the hyper-sphere § n is given by 

<->n — 



r(2±i) : 



where T is the well known gamma function, defined as 

T(x) = / t"- 1 e-*tft. 



The acceptance probability is 



Using the fact that [H Ex. 9.44] 



and in particular for any n > 1, 



p(ok) 



2 r(ffi) 

^ r(f) • 



r(2±i) fir _ i 



1 ^W) /n 

2V 2 - r(§) - V 2 
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Therefore 



The entropy of the messages is given by 



— - <p(ok) < J— . 

2im V 



log ( - J + - — - log 



1-p 



1 ., 

Inn 



H(P P ) < log(V27rn) + log 



' 1 1 / 2 



So in our case, we get 



27rra \ y 7rn 

< ^log(n) + 0(l). 

So, with n = 2c? 2 — 1, log(d) + 0(1) bits of communication are sufficient on average to simulate the joint quantum 
correlations of an arbitrary bipartite qudit pair, measured according to a TBO. □ 



IV. DISCUSSION AND CONCLUSION 



We have shown that in the general case of bipartite qudit pairs, we can apply the sampling approach of |9( to 
simulate the joint quantum correlations, using log(d) bits of communication on average. 

There are very few results in the literature concerning settings that are more general than projective measurements 
on maximally entangled qubit pairs. In the case of qudits, Bacon and Toner have shown how to simulate joint 
quantum correlations arising from TBO measurements with constant communication on av erag e, but the correlations 
are simulated approximately |l2j |. whereas here we simulate the correlations exactly. In |lOj |. Massar et al gave a 
protocol that simulates the correlations of any local measurement on an arbitrary bipartite state exactly, but within 
a different model, which uses communication only and no local hidden variables. 

In 0, we considered two ressources other than communication: post-selection, and nonlocal boxes. In the qubit 
setting, instead of iterating the rejection method until a suitable sample was selected, it could be proven that if the 
first sample failed the selection, then the second sample could be used. In this case, one bit of communication sufficed 
for Alice to communicate to Bob which sample to use, and the method could also be adapted to obtain a protocol 
that makes a single use of a nonlocal box. 

It turns out that this so-called "choice method" does not extend directly to dimensions other than two. Therefore, 
we do not obtain a worst case analysis in this more general setting, nor do we get a protocol using nonlocal boxes. 

On the other hand, our analysis immediately applies to protocols using post-selection, that is, where the protocol is 

allowed to abort with some probability. Here, the protocol succeeds with probability at least \J where n = 2d 2 — l, 
that is, <3(±). 
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